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FOREWORD 


As  part  of  the  on-going  research  program  in  "Decision  and  Control 
Models  in  Operations  Research,  Mr.  John  G.  Klincewicz  has  investigated 
the  behavior  of  multi-item  inventory  control  systems  in  which  the 
underlying  demand  distributions  are  estimated  from  a limited  sample  of 
historical  observations.  His  work,  which  carries  on  from  that  of 
Mr.  Richard  Ehrhardt  (Technical  Report  No.  7),  utilizes  demand  distribu- 
tions with  constant  standard  deviation-to-mean  ratio  of  unity.  Ehrhardt 's 
previous  experiments  with  Power  Approximation  control  considered  systems 
with  constant  variance- to-mean  ratios.  To  assist  the  reader  in  compre- 
hending the  results,  Klincewicz  presents  his  findings  in  tables  that 
also  exhibit  the  corresponding  results  of  Ehrhardt;  hence,  this  report 
is  nearly  self-contained. 


Other  related  reports  dealing  with  this  research  program  are  given 


below. 
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ABSTRACT  | 

THE  POWER  APPROXIMATION:  CONTROL  OF  MULTI-ITEM  INVENTORY  SYSTEMS 
WITH  CONSTANT  STANDARD-DEVIATION-TO-MEAN  RATIO  FOR  DEMAND 

John  G.  Kllncewlcz 

In  this  simulation  study,  we  consider  the  performance  of  the  Power 
Approximation  [Ehrhardt  (1976)]  for  a multi-item  inventory  system  in  which 
the  underlying  demand  distributions  are  negative  binomial  with  constant 
standard ^eviation-to-mean  ratio  of  1. 

The  Power  Approximation  formulas  [Ehrhardt  (1976)]  were  obtained  by 
least  squares  regression,  using  data  from  Inventory  Items  with  Poisson 
deznand  distributions  and  negative  binomial  demand  distributions  with 
variance-to-mean  ratios  of  3 and  9.  Previous  experiments  have  all  con- 
sidered the  effectiveness  of  the  Approximation  for  systems  with  constant 
varlance-to-mean  ratio. 

We  examine  the  constant  standard-deviation-to-mean  system  for  the 
situation  of  full  information  and  for  the  situation  in  which  the  decision-  ; 

maker's  knowledge  is  limited  to  a sample  of  previously-realized  demands.  | 

In  addition,  the  research  examines  the  accuracy  of  statistical  forecasts  ] 

that  predict  the  future  behavior  of  the  operating  characteristics.  As  a | 

result,  an  inventory  systems  designer  is  apprised  of  both  the  costs  of  1 

imperfect  information  and  the  extent  of  bias  in  the  forecast  estimates. 

In  general,  the  behavior  observed  for  the  constant  standard-deviation-to- 
mean  system  is  comparable  to  the  behavior  of  the  previously  studied  systems 


with  constant  varlance-to-mean  ratios  of  9. 
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1.  DESIGN  OF  INVENTORY  CONTROL  POLICIES 

1.1  The  Model 

Throughout  this  report,  we  deal  with  a single-item  Inventory  model. 

We  assume  periodic  review  of  an  Item's  Inventory  level  and  employ  a 
stationary,  discrete-time  stochastic  process  to  generate  the  Item's 
demand.  The  demand  sequence  consists  of  Independent 

Identically  distributed  negative  binomial  random  variables  taking  non- 
negative  Integer  values. 

Demands  are  met  as  long  as  stock  on  hand  Is  sufficient;  when  a 
stockout  occurs,  the  unfilled  demand  Is  completely  backlogged  until  a 
stock  replenishment  eventually  arrives. 

Items  kept  In  Inventory  are  assumed  to  be  conserved,  there  being 
no  losses  by  deterioration,  obsolescence,  or  pilferage;  disposal  Is  not 
allowed.  Inventory  on  hand  at  the  end  of  a current  period  Is  the  Inventory 
from  the  previous  period  plus  any  replenishment  that  arrives,  less  demand 
In  the  current  period.  Negative  Inventory  on  hand  represents  the  amount 
of  backlogged  demand.  Replenishments  are  assumed  to  be  delivered  a fixed 
lead  time  X periods  after  being  ordered.  The  time  sequence  of  events  In 
any  period  Is  taken  to  be  order,  delivery,  demand. 

We  assume  no  time  discounting  of  costs  and  postulate  an  unbounded 
horizon  over  which  the  Item  Is  demanded  and  stocked.  We  seek  Co  minimize 
expected  total  cost  per  period. 

The  cost  of  a replenishment  quantity  q Is  assumed  linear  with  fixed 


ordering  cost  K and  constant  unit  cost  c 


(1) 


c(q) 


K + cq 
0 


for  q > 0 
for  q « 0 , 


2 


Since  Items  are  not  lost  from  Inventory,  demand  Is  completely  filled,  and 
costs  are  not  discounted,  the  constant  unit  cost  c Is  not  a factor  In 
choosing  a minimum  cost  policy,  and  Is  suppressed  hereafter. 

The  Inventory  holding  cost  Is  proportional  to  any  stock  on  hand  at 
unit  cost  h 


1 

hi 

for  1 > 0 

(2) 

h(l)  = 

1 

0 

< 

for  1=0  ^ 

and 

the  unit  penalty  cost  ir 

Is  applied  to  any  quantity  on  backorder  at 

the 

end  of  each  period 

I ° 

for  1^0 

(3) 

tt(1)  = 

— it1 

for  1 < 0 ^ 

The  resulting 

total  cost 

function,  therefore.  Is  linear  In  K , tt  , 

and  hf  and  we  may  scale  these  parameters  so  that  the  value  of  the  unit 
holding  cost  h Is  unity.  Non-trlvlal  changes  In  cost  arise  only  with 
changes  In  the  ratios  K/h  and  ir/h  . 

We  postulate  that  control  over  replenishment  Is  exercised  by  an 
(s,S)  policy;  whenever  Inventory  x on  hand  and  on  order  at  the  start 
of  a period  drops  below  the  value  "little  s,"  an  order  Is  placed  for  a 
replenishment  of  size  S-x  . 

Given  our  assumptions,  when  the  demand  distribution  and  the  economic 
parameters  are  known,  there  Is  an  optimal  policy  that  has  the  (s,S)  form 


[Iglehart  (1963a, b),  Veinott  & Wagner  (1965)].  When  the  demand  distrlbu 


tion  is  not  known,  even  though  this  is  the  only  assumption  relaxed,  an 


optimal  policy  may  no  longer  be  of  the  (s,S)  form.  Nevertheless,  in 


) policy,  since  It  is  in  popular  use 


tnls  experiment  we  employ  an 


in  the  applied  situation  of  incomplete  information 


1.2.1  The  Power  Approximation 


The  Power  Approximation  devised  by  Ehrhardt  (1976)  is  an  algorithm 


for  computing  approximately  optimal  values  for  (s,S)  using  only  the  mean 
U and  variance  of  demand.  The  algorithm  is  executed  as  follows.  Let 


where  U(z)  is  given  by 


U(z)  = .182/z  + 1.142  - 3.466Z 


Otherwise 


compute 


4 


where  v is  the  solution  to 


♦ (v)  - ir/(ir+h)i 


is  the  cumulative  distribution  function  of  the  unit  normal 


distribution.  The  policy  parameters  are  then  given  by 


s “ minimum  {s 


S - minimum  {S 


If  demands  are  integer  valued,  s 


nearest  integer.  We  analyze  this  algorithm's  performance  in  Section  2 


Of  course,  in  real  applications,  the  mean  and  variance  of  demand 


are  not  always  known.  For  the  situation  where  only  sample  statistics  of 


previous  deiiiands  are  available,  the  Statistical  Power  Approximation 


described  below  can  be  implemented.  This  decision  rule  implicitly  assiimes 


that  demand  is  stationary.  Such  a rule  derived  for  stationary  conditions 


may  be  a reasonable  approximation  to  an  optimal  rule  when  the  demand 


process  is  mildly  nonstationary,  provided  that  the  policy  parameters 


and  S are  revised  periodically  to  meet  the  changing  conditions 


We  assume  in  this  study  that  a demand  history  of  fixed  length  is 


kept  to  make  each  revision,  and  equal  weight  is  given  to  each  observation 


This  is  not  optimal  if  the  demand  process  is  known  to  be  stationary,  for 


then  the  entire  history  should  be  accumulated  to  give  prngiessively  better 


knowledge  and  performance.  Even  when  demand  Is  known  to  be  nonstationary 


but  varying  In  a regular  manner,  such  as  by  a trend  or  periodic  cycle,  or 


both,  an  optimal  decision  rule  would  generally  utilize  the  entire  history 


5 

The  decision-maker  usually  Is  not  In  a position  to  know,  however,  that 
conditions  observed,  even  over  the  entire  past  history,  will  continue  to 
prevail.  This  provides  justification  for  making  frequent  revisions,  placing 
greater  weight  on  observations  from  the  Immediate  past  and  less  on  earlier 
history.  For  this  study  the  admittedly  arbitrary  choice  has  been  made  to 
keep  a history  of  fixed  length  and  give  equal  weight  to  all  observations  In 
this  history.  Let  T be  the  number  of  periods  between  policy  revisions, 
which  will  be  termed  the  revision  Interval;  assume  that  a history  of  T 
periods'  demands  Is  kept  for  use  at  each  revision. 

The  statistics  required  by  our  decision  rules  are  the  sample  mean  and 
variance  of  demand,  5 and  v , respectively.  If  t Is  a period  at  the 
! beginning  of  which  revision  Is  made,  then 


I 

I 


(7) 


^ 111  ^t-r 

V . (T-l)-l  . 


When  using  the  Statistical  Power  Approximation  we  periodically  obtain 
values  for  (s,S)  by  substituting  C and  v for  y and  In  equations 

(1)  through  (6).  We  Investigate  the  properties  of  the  Statistical  Power 
Approximation  In  Sections  3 and  4 of  this  study. 

1.2.3  Design  Parameters 

In  the  basic  experiment,  each  run  Is  made  with  a different  combination 
of  parameters.  A negative  binomial  distribution  Is  used  to  generate  the 
demand  sequences  with  standard-devlatlon-to-mean  ratio  fixed  at  1.  Runs 


« i 


are  made  with  four  mean  values  for  the  distributions,  shown  In  Table  1.1. 


r 


Three  values,  0,  2,  and  4,  are  assigned  to  lead time.  Since  the  cost 
function  Is  linear  In  the  parameters  K , h , and  tt  , the  value  of  the 
unit  holding  cost  h Is  set  at  unity.  Essential  changes  In  cost  arise 
only  with  changes  In  the  ratios  K/h  and  ir/h  . The  stockout  costs  are 
TT  * 4,  9,  and  99  and  the  fixed  cost  values  are  K = 32  and  64  . All 
combinations  of  Input  levels  are  Included  In  the  system  corresponding  to 
a full  factorial  design.  Table  1.1  summarizes  the  parameter  specifications 
of  the  system. 

The  resulting  combinations  of  values  for  mean,  leadtime,  stockout 
cost,  and  fixed  ordering  cost  total  72.  Each  set  of  Input  parameters  is 
run  with  a revision  Interval  and  revision  history  length  of  26  periods,  using 
demand  data  accumulated  over  the  same  time  Interval.  In  each  run  of  the 
experiment,  200  revisions  and  forecasts  are  made. 


Table  1.1 
System  Parameters 


Factor 


Levels 


Demand  distribution 

Mean  demand 

Unit  holding  cost 

Unit  backlog  penalty  cost 

Replenishment  setup  cost 

Replenishment  leadtime 


Negative  Binomial  (a/p 
2,  4,  8,  16 
1 

4,  9,  99 
32,  64 
0,  2,  4 


= 1) 


Number  of  Levels 


3 

2 

3 


For  a detailed  discussion  of  demand  generation,  system  Initialization,  and 
output  analysis,  see  MacCormick  [(1974,  pp.  52-58)]. 


1.2,4  Item  Operating  Characteristics 


Detailed  operating  characteristics  of  each  Item  in  the  system  are 
accumulated.  These  are  the  period-end  inventory  on  hand,  the  period-end 
stockout  quantity,  the  frequency  of  period-end  stockouts,  the  replenish- 
ment quantity,  the  frequency  of  replenishment,  and  the  total  cost  incurred 
It  is  assumed  that  the  decision-maker  is  interested  only  in  the  average 
performance  between  revisions,  corresponding  to  the  practice  adopted  by 
accountants  of  making  reports  to  management  at  periodic  intervals. 

One  objective  of  the  experiment  is  to  make  inferences  about  the 
distributions  of  these  operating  characteristics  for  each  of  the  three 
decision  rules  (optimal  policies  under  full  information,  approximately 
optimal  policies,  and  statistical  approximation  policies).  Exact  expected 
values  are  computed  analytically  for  each  operating  characteristic  under 
optimal  control  using  full  information  and  under  approxl'iiately  optimal 
control  using  full  information.  Under  statistical  control  the  simulation 
experiment  yields  a sequence  of  200  sample  values  of  each  operating  charac 
terlstic  for  which  confidence  intervals  are  constructed  using  the  methods 
discussed  in  MacCormlck  (1974) . 

1.2.5  Forecasting 

The  forecasting  method  used  here  is  the  same  as  that  described  by 
MacCormlck  (1974) , namely,  retrospective  simulation.  This  method  takes 
the  recent  demand  history  for  the  items  in  the  sample  and  estimates  the 
performance  of  the  chosen  policy  for  this  history.  Thus,  the  method  uses 
the  same  data  twice,  once  to  fix  the  policy  parameters,  and  once  to  fore- 
cast the  performance.  As  a result,  the  forecasts  are  biased,  and  this 
study  investigates  the  extent  of  the  bias. 
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At  each  revision,  after  the  (s,S)  parameters  have  been  set  to  new 
values,  a forecast  Is  made  of  the  properties  of  the  system.  A history  of 
T+X  demand  values  Is  kept  to  make  the  forecast.  [We  note  that  there  are 
differences  In  Information  storage  requirements  of  considerable  Importance 
In  practice.  To  make  a forecast  of  performance  by  retrospective  simulation, 
the  demands  for  the  Item  in  each  of  T+X  periods  must  be  stored,  whereas 
the  statistical  decision  rule  may  require  only  a handful  of  sufficient 
statistics  to  set  (s,S).  Also,  for  multi-item  systems  the  cost  of  accumu- 
lating in  storage  the  data  needed  to  make  forecasts  by  retrospective 
simulation  is  likely  to  be  so  high  that  histories  will  be  kept  for  only  a 
representative  sample  of  items  from  the  system.] 

The  forecasts  are  made  by  running  the  system,  using  the  new  values  of 
s and  S , and  observing  the  operating  characteristics  as  if  these  (s,S) 
values  had  been  in  force  when  the  history  occurred.  Each  time  there  is  a 
policy  revision,  there  is  a forecast  for  each  operating  characteristic. 

The  simulation  initializes  the  stock  on  hand  and  on  order  at  the  actual 
value  at  the  time  of  the  forecast,  that  is,  at  the  time  of  revision.  As 
actual  stock  on  hand  is  not  recorded  until  the  elapse  of  a leadtime  after 
revision,  the  same  interval  is  allowed  to  elapse  before  recording  it  for 
the  forecast.  The  inventory  on  hand  variable  is  therefore  initialized  at 
the  initial  value  for  inventory  on  hand  and  on  order  less  the  first  X 
demands  in  the  history  of  T+X  demands. 

1.3  A Multi-Item  System 

Scientific  techniques  for  inventory  control  are  generally  applied  to 
systems  of  many  items.  This  experiment  combines  the  results  from  the  72 
single-item  simulations  into  a multi-item  system.  Since  management  gener- 
ally assesses  the  performance  of  control  techniques  for  a multi-item 


inventory  system  by  observing  Indices  that  are  aggregate  operating  charac- 
teristics [Wagner  (1962)],  certain  aggregate  characteristics  have  been 
computed. 

The  operating  characteristics  of  the  multi-item  system  under  statistical 
control  have  been  measured  by  aggregating  the  sample  values  of  the  corres- 
ponding characteristics  for  each  item  in  the  system.  When  the  system  is 
operated  under  perfect  information,  these  characteristics  are  computed 
analytically.  The  aggregate  of  average  total  cost  per  period  is  computed 
as  the  arithmetic  sum  of  the  corresponding  costs  for  each  item.  The 
components  of  total  cost  for  inventory  storage,  backlog  penalty,  and  for 
ordering  replenishments  are  similarly  computed.  The  aggregate  backlog  and 
replenishment  frequencies  are  arithmetic  averages  of  the  corresponding 
frequencies  observed  for  each  item  in  the  system.  Since  the  unit  inventory 
holding  cost  for  all  items  is  unity,  the  average  number  of  units  in  inventory 
at  period-end  is  identical  to  the  aggregate  average  holding  cost  per  period. 
Finally,  a weighted  proportion  of  demand  backlogged  is  computed  as  the  ratio 
of  a weighted  sum  of  the  average  quantity  backlogged  per  period  to  a 
weighted  sum  of  the  (exact)  mean  values  of  demand.  The  weights  used  in 
both  the  numerator  and  denominator  of  the  ratio  are  the  unit  cost  of 
backlogging  demand  for  the  respective  item. 

1 . 4 Use  of  Constant  Standard'^eviation-to-Mean  Ratio 

The  Power  Approximation  formulas  [Ehrhardt  (1976)]  were  obtained  by 
least  squares  regression,  using  values  of  s,S,  and  D = S-s  for  288  items 
as  data  points.  Three  types  of  demand  distributions  were  Included  in 
these  items:  Poisson  and  negative  binomial  with  variance- to-mean  ratios  of 
3 and  9.  Previous  experiments  have  all  considered  the  effectiveness  of  the 


s 


Table  1.3  displays  the  standard  deviation- to-mean  ratio  for  items 
with  mean  y » 2,  4,  8 and  16  when  the  variance-to-mean  ratio  equals 
3 and  9.  The  average  standard-deviation- to-mean  ratio  for  these  items 
is  1.07. 


Table  1.3 


Standard-Devlatlon-to-Mean  Ratios  for  Items  in 


Fixed  Variance-to-Mean  Systems 


Mean  y 


Standard  Deviation/Mean^g/y 


Variance/Mean  = 9: 


Variance/Mean  = 3: 


Average  = 1.36 


Average  = .78 


„ ...  .78  + 1.36  , 

Overall  Average  = r = 1.07 
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2.  SYSTEM  CONTROL  WITH  FULL  INFORMATION  ABOUT  DEMAND 

In  this  section  we  describe  the  Inventory  system  with  constant 
standard-devlation-to-mean  ratio  under  optimal  and  Power  Approximation 
control  with  full  information  about  demand.  We  seek  to  establish  bench- 
mark values  of  the  system  operating  characteristics  for  later  comparison 
with  a system  controlled  using  only  statistical  information  about  demand. 

2 . 1 Optimal  Control  with  Full  Information 

Each  item  in  the  72  item  system  described  in  Section  1 is  controlled 
using  an  optimal  (s,S)  policy.  Optimal  values  for  the  policy  parameters 
are  computed  using  the  algorithm  in  Velnott  and  Wagner  (1965) . The  expected 
values  of  the  operating  characteristics  are  calculated  by  first  determining 
their  conditional  expectations  for  given  values  of  inventory  on  hand  plus 
on  order  after  ordering.  The  distribution  of  inventory  on  hand  plus  on 
order  after  ordering  is  then  used  to  form  the  unconditional  expectations. 

Table  2.1  lists  the  resulting  expected  values  of  average  total  cost  per 
period  and  its  components  for  the  72  item  negative  binomial  system  with 
standard-deviatlon-to-mean  ratio  of  1 and  the  72  item  systems  with 
variance-to-mean  ratios  of  3 and  9.  Each  component's  percent  of  total 
cost  is  shown  in  parentheses.  For  all  quantities  shown  in  the  table,  the 
value  associated  with  the  a/\i  = 1 system  is  between  the  comparable 
values  for  the  a^/p  = 9 and  o^/p  = 3 systems.  In  general,  however, 
the  values  are  closer  to  those  for  the  higher  variance-to-mean  ratio 
system. 
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Table  2.1 

Average  Costs  per  Period  for  a 
Multi-Item  System  Optimally  Controlled 


Cost  Component 

Standard 

Deviation/Mean=l 

Variance/Mean=9 

Varlance/Mean=3 

Inventory 

1897  (59.9) 

1935  (59.7) 

1273  (54.4) 

Backlog 

629  (19.9) 

682  (21.0) 

343  (14.6) 

Replenishment 

642  (20.3) 

626  (19.3) 

728  (31.1) 

Total 

3169  (100.0) 

3243  (100.0) 

2345  (100.0) 

The  apportionment  of  expected  average  costs  for  various  classifications 
of  the  items  in  each  multi-item  system  is  shown  in  Table  2.2,  as  percentage 
distributions  of  the  estimated  expected  average  total  cost  per  period.  We 
note  that  items  with  mean  = 16  contribute  a much  greater  proportion  of  total 
costs  in  the  o/y  = 1 system  than  in  either  the  o^/y  =9  or  the  a^/y  = 3 
systems.  This  is  due  to  the  fact  that  items  with  y * 16  have  a higher 
variance  for  the  fixed  standard-deviatlon-to-mean  system  than  in  either  of 
the  fixed  varlance-to-mean  systems  (see  Section  1.4). 

Aggregate  values  of  the  other  operating  characteristics  for  the  various 
classifications  are  set  out  in  Table  2.3.  The  classifications  are  one-way 
in  the  sense  that  items  in  the  system  are  grouped  according  to  the  value 
taken  by  a single  input  parameter. 

The  sensitivity  of  expected  costs  per  period  and  other  operating 
characteristics  is  summarized  in  Table  2.4,  which  has  been  constructed  by 
observing  any  monotonlclty  in  the  relation  for  the  multi-item  system  (see 
Table  2.2).  We  note  two  differences  between  the  patterns  observed  for  the 
a/y  “ 1 system  and  those  for  the  o^/y  ■ 3 and  9 systems.  For  the 


Apportionment  of  Aggregate  Costs  per  Period  for  a 72-Item  System  Under 
Optimal  Control  with  Full  Information 
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Operating  Characteristics  of  a 72-Item  System  Under 
Optimal  Control  with  Full  Information 


REPLENISHMENT  FREQUENCY  1.222  .213  .222  .230  .254  .189  .235  .220  .210 
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a^/y  ■ 9 system,  proportion  demand  backlogged  decreases  monotonically 
with  respect  to  setup  cost,  whereas  for  the  a^/y  “ 3 system  the 
relationship  is  monotonically  Increasing.  For  the  a/y  = 1 system  the 
relationship  Is  essentially  constant. 

There  is  no  discernible  pattern  in  the  relationship  of  expected 
backlog  cost  with  respect  to  penalty  cost  for  the  constant  standard- 
devlatlon-to-mean  system.  There  is,  however,  a monotonically  increasing 
relationship  in  the  variance/mean  = 9 system  and  a monotonically  decreasing 
relationship  for  the  variance/mean  = 3 system. 

The  previous  experiments  of  MacCormick  (1974)  and  Estey  and  Kaufman 
(1975)  for  negative  binomial  systems,  and  of  Kaufman  and  Klincewicz  (1976) 
for  a sporadic  demand  system,  with  constant  variance-to  mean  ratio  have 
Indicated  that  the  dependence  of  expected  average  total  cost  per  period 
on  the  input  parameters  appears  to  be  multlplicatively  separable  in  form. 

Table  2.4  shows  the  percentage  distributions  according  to  input 
parameter  groupings  of  expected  total  cost  per  period  for  the  72  item 
system  with  a/y  ” 1 . There  are  three  entries  in  each  cell  of  Table  2.4. 
The  first  is  the  proportion  of  expected  total  cost  per  period.  Values  in 
the  second  entry  are  calculated  as  the  product  of  the  corresponding 
proportions  observed  for  the  one-way  distributions  that  form  the  border 
to  Table  2.5.  For  example,  when  mean  = 2 and  ir/h  ■ 99,  the  product 
is  .096 (.459)  = .044  , which  compares  to  the  observed  proportion  of  .040 
in  the  two-way  distribution.  The  third  entry  in  the  cell  is  the  difference 
between  the  observed  two-way  proportion  and  the  product  of  the  observed 


one-way  proportions. 
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Table  2.4 


Sensitivity  to  Input  Parameters  of  Operating  Characteristics 
of  an  Optimally  Controlled  Inventory  System 
with  Negative  Binomial  Demanas 


INPUT  PARAMETERS 

Mean 

K/h  Lead  time  Demand 


STANDARD  DEVIATION/MEAN 


OPERATING  CHARACTERISTIC 


E(COST) 


t i + 


+ + i 


+ + + + i + 


E (PERIOD-END  INVENTORY) 


+ + + 


11+  t + + 


E (BACKLOG  COST) 


? i + 


+ + + 


+ + i + + + 


E (BACKLOG  FREQUENCY) 


+ + -I- 


+ - + 


+ + + + i i 


PROPORTION  DEMAND  BACKLOGGED 


+ + + 


- 4-  + 


+ + + +4-1 


E (REPLENISHMENT  COST) 


+ + + 


+ 4 4 


4 4 4 4 4 4 


E (REPLENISHMENT  FREQUENCY) 


4 4- 


4 4 4 


4 4 4 4 4 4 


POLICY  PARAMETERS 


444  444  444 


4 4 4 


44  444  444 


4 4 4 


444  444  444 


The  three  symbols  in  each  column  refer  to  the  following  three  sets  of  demand 


distributions:  Negative  Binomial  (standard  deviation/mean  = 1),  Negative 


Binomial  (Variance/Mean  = 9)  and  Negative  Binomial  (Variance/Mean  = 3) 


respectively. 


4 indicates  a monotonically  increasing  relationship 
4 indicates  a monotonically  decreasing  relation«;h ip 
- indicates  an  approximately  constant  relationship 
? indicates  that  no  simple  relationship  could  be  discerned 
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The  observed  differences  are  all  relatively  small  and  give  evidence 
of  a multiplicative  form  for  the  dependence  of  expected  total  cost  per 
period  on  input  parameters  for  the  constant  o/y  case  as  well.  As  in  the 
fixed  varlance-to-mean  ratio  case,  we  observe  trends  in  the  differences 
with  respect  to  changes  in  certain  parameter  settings.  The  trends  that  are 
observed  across  the  rows  for  mean  demand  values,  however,  are  all  opposite 
to  the  trends  observed  for  the  constant  variance-to-mean  case.  For  instance, 
in  Table  2.5  we  note  that  for  classifications  with  y = 2 , the  differences 
change  by  negative  Increments  as  we  consider  higher  values  of  penalty  cost 
IT  . For  constant  variance-to-mean  ratios,  these  increments  were  positive. 
Further  experiments  may  show  if  these  observed  trends  extrapolate  beyond 
the  range  of  parametars  considered  here. 

Finally,  we  explore  an  approximation  to  the  backlog  frequency  that  is 
well  known  to  writers  on  inventory  control.  For  an  optimally  controlled 
multi-period  Inventory  model  with  linear  costs,  the  ratio  l/(l+n)  forms 
an  upper  bound  to  the  probability  of  backlog.  Table  2.6  shows  the  backlog 
frequencies  are  at  or  below  the  upper  bounds  for  all  three  inventory  systems 
we  are  considering. 


Table  2.6 

Backlog  Frequencies  for  a Multi  Item  System 
with  Negative  Binomial  Demands 


4 

TT/h 

9 

99 

l/(l+7r) 

.2 

.1 

.01  ? 

i 

Standard  Deviation/Mean  “ 1 

.186 

.094 

.009  - 

Variance/Mean  • 9 

.189 

.095 

.010 

Variance/Mean  - 3 

.182 

.092 

.009 

J 
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2 . 2 Power  Approximation  Control  with  Full  Informati-jn 

We  obtain  a Power  Approximation  policy  for  each  Item  In  our  multi-item 
system  using  equations  (1)  through  (6)  In  Section  1.2.3.  Although  the  control 
Is  described  as  utilizing  full  Information,  only  the  mean  and  variance  of  the 
demand  distribution  need  to  be  known  In  order  to  compute  (s,S). 

Table  2.7  lists  the  resulting  expected  value  of  average  total  cost 
per  period  and  Its  components  for  the  72  Item  system  with  a/p  = 1 , with 
a^/p  =9,  and  with  a^/p  = 3 . Each  component's  percent  of  total  cost 
Is  shown  In  parentheses. 


Table  2.7 

Average  Costs  per  Period  for  a Multi-Item  Negative 
Binomial  System  with  Power  Approximation  Control 


COST  COMPONENT 

Standard 

Deviation /Mean  = 1 

Variance/Mean  = 9 

Variance/Mean  = 

INVENTORY 

1939  (61.0) 

1962  (60.3) 

1262  (53.7) 

BACKLOG 

591  (18.6) 

664  (20.4) 

362  (15.4) 

REPLENISHMENT 

651  (20.5) 

626  (19.3) 

726  (30.9) 

TOTAL 

3181  (100.0) 

3252  (100.0) 

2350  (100.0) 

The  values  for  each  component  of  total  cost  for  the  standard  deviation/ 
mean  > 1 case  fall  between  the  values  for  the  variance/mean  ° 9 and  the 
variance/mean  ••  3 cases.  As  Is  the  case  under  optimal  control,  however,  the 
results  are  much  closer  to  the  variance/mean  = 9 results,  which  show  greater 
total  cost  than  the  variance/mean  • 3 results. 
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In  Table  2.8  we  list  absolute  and  percentage  differences  between  optimal 
and  Power  Approximation  control.  We  see  that  the  fixed  standard  devlatlon/mean 


system  shows  only  slightly  worse  degradation  In  total  cost  than  the  fixed 
variance/mean  cases.  The  percentage  Increase  In  Inventory  holding  costs 
and  replenishment  cost  and  the  percentage  decrease  In  backlog  cost, 
compared  to  optimal  control, are  more  pronounced  for  the  o/y  = 1 case 
than  for  the  a^/y  = 9 case.  For  the  a^/y  = 3 case.  It  Is  the  backlog 
costs  which  are  above  those  for  optimal  control  whereas  the  Inventory  and 
replenishment  costs  show  percentage  decreases. 

Table  2.8 

Average  Costs  Per  Period  for  a Multi  Item  Negative  Binomial 
System;  Comparison  of  Optimal  and  Power  Approximation  Control 


COST  COMPONENT 

Standard 

Devlatlon/Mean  ® 1 

Variance/Mean  = 9 

INVENTORY 

42  (2.2) 

26  (1.4) 

BACKLOG 

-39  (-6.1) 

-18  (-2.7) 

REPLENISHMENT 

9 (1.4) 

1 (0.2) 

TOTAL 

12  (0.4) 

9 (0.3) 

-13  (-1.0) 
17  (5.0) 
-1  (-0.1) 


3 (0.1) 


Note:  Table  2.8  shows  absolute  Increase  or  decrease  In  the  cost  components  of 
Table  2.7  over  those  In  Table  2.1  and  percentage  changes  In  parentheses. 


The  apportionment  of  aggregate  costs  per  period  for  various  parameter 
classifications  Is  shown  In  Table  2.9,  which  can  be  compared  with  Table  2.2. 
The  differences  in  apportionment  of  aggregate  costs  between  the  constant 
standard  devlatlon/mean  system  and  the  constant  variance/mean  systems, 
i^lch  can  be  seen  In  Table  2.9,  reflect  the  differences  In  apportionment 
under  optimal  control  as  seen  In  Table  2.2.  Detailed  comparison  of  the 
three  systems  under  optimal  and  approximately  optimal  control,  expanding 
on  Table  2.8,  Is  In  Table  2.10  showing  percentage  Increases  In  each  component 


Apportionment  of  Aggregate  Costs  per  Period  for  a 72-Item  System  Under 
Power  Approximation  Control  with  Full  Information 


Table  2.10 
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of  cost . 

The  percentage  decrease  In  backlog  cost  and  percentage  Increase  In 
Inventory  cost  from  the  optimal  that  were  noted  In  Table  2.8  fur  the 
o/y  1 system  are  most  noticeable  for  the  subsystems  with  parameter 
values  TT  = 99  and  y = 16  . The  change  In  total  cost  from  the  optimal, 
however.  Is  not  greater  than  0.6%  for  any  subsystem.  Overall,  the  Power 
Approximation  with  full  Information  appears  to  perform  close  to  optimal 
for  the  constant  standard-devlatlon-to-mean  system. 

Table  2.11  shows  the  values  of  other  operating  characteristics  for 
the  systems  under  approximately  optimal  control.  Percentage  changes  In 
the  values  of  the  operating  characteristics  In  Table  2.11  from  those  In 
Table  2.3  are  set  out  In  Table  2.12.  The  percentages  from  optimality  for 
proportion  demand  backlogged  with  parameter  settings  ir  = 99  and  y = 16 
reflect  the  pattern  noted  earlier  for  backlog  cost.  Overall,  for  the 
o/y  • 1 system,  the  period-end  Inventory  and  replenishment  frequency 
Increase  and  the  proportion  demand  backlogged  decreases  from  the  optimal 
values  as  they  do  for  the  a^/y  * 9 system;  backlog  frequency  Increases 
above  optimal  overall,  as  It  does  for  the  a^/y  = 3 system. 

Table  2.13  summarizes  the  sensitivity  to  the  values  of  input 
parameters  shown  by  the  operating  characteristics  of  a negative  binomial 
system  with  standard-devlati jn-to-uean  of  1 and  with  variance-to-mean  of 
9 and  of  3 under  Power  Approximation  control.  No  major  discrepancies 
between  the  pattern  for  the  constant  standard  deviation/mean  case  and  the 
pattern  for  the  constant  variance/mean  case  are  evident.  We  note,  however, 
the  sensitivity  of  the  backlog  characteristics  with  respect  to  leadtime  and 
of  the  backlog  cost  with  respect  to  penalty  cost.  For  these  relationships, 
no  clear  pattern  Is  discernible  for  o^/y  ■ 9;  the  patterns  observed  for 


Percentage  Excess  of  Operating  Characteristics  for  a 72-Item  System  Under  Power 
Approximation  Control  over  the  Same  Characteristics  for  the  System  Under  Optimal  Control 
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Table  2.13 

Operating  Characteristics  of  a 72-Item  Inventory  System 
Using  the  Power  Approximation:  Sensitivity  to  Input  Parameters 


INPUT  PARAMETERS 

Mean 

K/h  Leadtime  Demand 


OPERATING  CHARACTERISTICS 
E(COST) 

E (PERIOD-END  INVENTORY) 

E (BACKLOG  COST) 

E (BACKLOG  FREQUENCY) 

PROPORTION  OF  DEMAND  BACKLOGGED 
E (REPLENISHMENT  COST) 

E (REPLENISHMENT  FREQUENCY) 


+ + f + + + + + + + + + 

+ +++  f + + 

+ ?+  + + + +?+  + 

+ + 4-  t + 4 -?-  ??- 

4 + 4 + + + +?+  444 

4 + 4 444  444 

444  444  444 


— — — 444  444  444 

444  444  444  444 

444  444  444  444 

The  three  symbols  in  each  column  refer  to  the  following  three  sets  of  demand 

distributions: (Standard  Deviation/Mean  = 1),  Negative  Binomial  (Variance/ 

Mean  - 9),  and  Negative  Binomial  (Variance/Mean  = 3),  respectively,  and 

4 indicates  a monotonlcally  increasing  relationship; 

4 indicates  a monotonlcally  decreasing  relationship; 

- indicates  an  approximately  constant  relationship; 

? indicates  that  no  simple  relationship  could  be  discerned. 


POLICY  PARAMETERS 
D 
s 
S 


o/u  « 1 in  these  cases  are  Identical  to  those  observed  for  o^/y  = 3. 

Finally,  we  examine  the  ratio  1/(1+tt)  as  an  approximation  to  the 
probability  of  backlog  in  any  period  for  each  item.  The  entries  in 
Table  2.1A  have  been  extracted  from  Table  2.11,  The  values  for  the 
approximation  are  quite  close  to  actual  for  the  system  under  study,  and 
quite  comparable  to  the  results  for  the  constant  variance-to-mean  ratio. 

Table  2.14 

Backlog  Frequencies  for  a Multi-Item  System 
with  Negative  Binomial  Demands 

Penalty  coat  ■ ir  /h 


4 

9 

99 

l/d-Hr) 

.2 

.1 

.01 

Standard  Deviation/Mean  = 1 

.200 

.091 

.008 

Variance/Mean  = 9 

.189 

.088 

.010 

Variance/Mean  = 3 

.204 

.097 

.009 
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3.  SYSTEM  CONTROL  WITH  STATISTICAL  DEMAND  INFORMATION 

We  suppose  that  a manager  responsible  for  the  control  of  the  Inventory 
system  discussed  In  this  section  assumes  that  demand  can  be  modeled  suffi- 
ciently well  by  a sequence  of  Independent  demands,  distributed  statlonarlly, 
or  approximately  statlonarlly.  But  all  that  the  manager  actually  knows 
about  demand  Is  the  Information  provided  by  a sample  from  the  realized 
demand  history.  The  sample  can  be  used  to  estimate  central  moments  of  the 
demand  distribution;  such  estimates  take  the  place  of  the  true  distribution 
moments  In  control  policies  derived  to  be  approximately  optimal  for  the 
case  of  perfect  Information.  Details  of  the  control  policies  were  given 
previously  In  Section  1. 

The  operating  characteristics  describing  the  performance  of  Inventory 
systems  under  control  using  only  statistical  Information  about  demand  have 
not  been  computed  analytically.  The  results  In  this  section  are  estimates 
of  the  performance  of  some  selected  Inventory  systems  found  by  computer 
simulation.  During  each  simulation  run  the  control  policy  parameters  are 
revised  at  regular  Intervals  to  allow  for  potential  changes  In  the  under- 
lying demand  distribution.  Such  changes  do  not,  In  fact,  occur  during 
the  experiments,  but  In  a real-life  system  periodic  revision  Is  essential 
to  detect  and  adapt  to  changing  demand  structure. 

We  compare  our  simulation  results  for  standard'nlevlatlon-to-mean 
ratio  of  1 with  those  of  Ehrhardt  (1976)  for  the  varlance-to-mean  ratios 
of  3 and  9.  Estimated  average  total  cost  per  period,  and  Its  components, 
are  shown  In  Table  3.1  for  each  of  these  three  systems. 

Table  3.2  shows  estimates  of  absolute  and  percentage  Incremental 
changes  In  expected  costs  per  period  when  statistical  control  Is  compared 


'fdaf 


'•‘^Kf.HC^.r*- 


30 


Table  3.1 

Average  Costs  per  Period  for  a 72-Item  System: 
Control  with  Statistical  Information  about  Demand 


Cost  Component 

Standard 

Deviation/Mean  = 1 

Variance/Mean  = 9 

Var iance/Mean  = 3 

Inventory 

Backlog 

Replenishment 

1990  (56.9) 
854  (24.4) 
652  (18.7) 

1989  (55.0) 
993  (27.5) 
634  (17.5) 

1297  (52.6) 
439  (17.8) 
732  (29.6) 

Total 

3496  (100.0) 

3616  (100.0) 

2467  (100.0) 

Table  3.2 

Average  Costs  per  Period  for 
Comparison  of  Statistical  Control 

a Multi  .Item  System: 
with  Optimal  Control 

Cost  Component 

Standard 

Deviation/Mean  = 1 

Variance/Mean  = 9 

Variance/Mean  = 3 

Inventory 

Backlog 

Replenishment 

91  (4.9) 
225  (35.8) 
9 (1.5) 

54  (2.8) 
311  (45.7) 
8 (1.3) 

22  (1.7) 
94  (27.4) 
4 (0.6) 

Total 

327  (10.3) 

373  (11.5) 

120  (5.1) 

Note:  Table  3.2  shows  the  absolute  increase  in  cost  components  of  Table  3.1 
over  those  of  Table  2.1  and  percentage  increases  shown  in  parentheses. 


with  optimal  control  given  full  information.  The  percentage  excess  in  total 
cost  for  the  standard-deviation-to-mean  of  1 system  (10.3%)  is  greater  than 
that  for  the  variance-to-mean  of  3 system  (5.1%)  but  slightly  less  than  the 
percentage  excess  for  the  variance-to-mean  of  9 system  (11.5%). 
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When  the  Items  In  each  system  are  dlvl'led  Into  subsystems  according 
to  the  Input  parameter  values,  we  obtain  the  apportionment  of  costs  among 
the  subsystems  shown  In  Table  3.3.  This  table  Is  similar  In  format  to 

Table  2,2,  which  shows  systems  controlled  with  full  information. 

Comparison  of  Table  3.3  and  Table  2.2  Indicates  that  a llgher  propor- 
tion of  expected  total  cost  Is  formed  by  the  backlog  component  when  the 
o/ii  = 1 system  Is  controlled  statistically.  This  same  pattern  was 
observed  for  the  constant  variance- to-mean  systems  as  v’ell. 

Table  3.4  shows  the  estimated  percentage  Incremental  changes  In  the 
expected  average  total  cost  and  Its  components.  For  the  various  subsystems 
this  table  elaborates  on  the  changes  summarized  In  Table  3.2.  In  comparison 
with  the  constant  varlance-to-mean  case,  the  y = 16  subsystem  shows  a 
greater  percentage  excess  In  total  cost  for  a/u  = 1 . As  noted  In 
Section  1.4,  the  variance  for  y = 16  with  o/y  =1  Is  greater  than  the 
variance  for  y = 16  with  either  o^/y  ■=  3 or  a^/y  = 9 . Likewise,  the 
u = 2 subsystem,  which  has  a smaller  variance  when  o/y  = 1 , shows  a 
smaller  percentage  excess  In  total  cost  when  o/y  = 1 than  when 
o^/y  * 3 or  9 . 

We  tabulate  In  Table  3.5  the  physical  operating  characteristics 
for  the  three  systems,  and  In  Table  3.6  we  set  out  the  percentage  Increase 
of  each  over  the  corresponding  values  In  Table  2.3.  There  Is  a general 
pattern  of  greater  percentage  Increase  In  period-end  Inventory  for  the 
constant  standard-devlatlon-to-mean  system  than  for  either  of  the  constant 
varlance-to-mean  systems. 

The  largest  percentage  excesses  for  the  o/y  > 1 case  occur  for 
backlog  frequency  (60. 4Z)  and  proportion  demand  backlogged  (78.7%)  In  the 
IT  ■ 99  subsystem.  These  values  are,  however,  not  as  large  as  for  the 


J 


Apportionment  of  Costs  for  a 72-Item  System  Controlled  with  Statistical  Information  about  Demand 
(Revision  Interval  26  Periods,  Revision  History  Length  26  Periods) 


Percentage  Excess  of  Costs  of  a 72-Item  System  under  Statistical  Control 
er  the  Costs  of  Controlling  the  Same  System  Optimally  with  Full  Information 
(Revision  Interval  26  Periods,  Revision  History  Length  26  Periods) 


Operating  Characteristics  of  a 72-Item  System  Controlled  with  Statistical  Information  About  Demand 
(Revision  Interval  26  Periods,  B^^'lslon  History  Length  26  Periods) 


DEMAND  BACKLOGGED  .022  .224  .090  .007  .020  .023  .015  .023  .028  .057  .036  .022  .014 

REPLENISHMENT  EHEQUENCY  .222  .223  .223  .221  .254  .191  .234  .219  .213  .119  .174  .248  .348 


Percentage  Excess  of  Operating  Characteristics  of  a 72-Item  System  Under  Statistical  Control 
Over  Those  for  System  Under  Optimal  Control  with  Full  Information 
(Revision  Interval  26  Periods,  Revision  History  Length  26  Periods) 


I 


I 


\ 
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a^/p  « 9 system  where  the  percentage  Increments  for  these  characteristics 
exceeded  100%. 

In  Table  3.7  we  Indicate  the  direction  of  sensitivity  of  the  operating 
characteristics  and  policy  variables  to  changes  In  Input  parameters  for 
each  of  the  three  systems  we  have  been  comparing.  The  only  discrepancy 
of  note  Is  the  fact  that,  for  the  o/y  * 1 system,  proportion  demand 
backlogged  shows  no  monotonlc  pattern  with  respect  to  mean  demand  whereas 
In  both  of  the  systems  with  a^/y  constant,  a monotonlc  decreasing  relation- 
ship Is  discernible. 

As  we  did  for  the  case  of  optimal  control  with  full  Information,  we 
consider  the  structural  form  of  the  dependence  of  expected  total  cost  per 
period  for  the  multl~ltem  system  on  the  levels  of  the  Input  parameters.  For 
optimal  control  in  Section  2 (Table  2.5)  It  was  observed  that  the  distribu- 
tion of  expected  total  costs  according  to  Input  parameter  values  closely 
conformed  to  a multiplicative  relationship.  Table  3.8  demonstrates  the 
fit  for  the  case  of  control  with  statistical  Information  about  demand, 
with  a 26  period  revision  Interval  and  utilizing  a 26  period  demand  history. 
The  differences  In  Table  3.8  are  closely  comparable  to  the  differences 
noted  In  Table  2.3.  The  same  trends  can  be  observed  by  scanning  the 
differences,  and  we  caution  against  extrapolation.  For  the  range  of 
parameters  tested,  however,  the  use  of  constant  standard-devlatlon-to- 
mean  ratio  has  not  Impaired  the  fit  of  the  multiplicative  structure  that 
was  first  observed  for  the  case  of  constant  varlance-to-mean  ratio. 

We  also  consider  the  ratio  1/(1+7t)  for  Its  value  as  an  approxima- 
tion to  the  frequency  of  backlog  for  a typical  Item  In  a multi- Item  system. 

As  with  the  constant  varlance-to-mean  systems,  the  actual  value  exceeds 
the  estimate  l/(l+ir)  for  all  three  values  of  tt  tested  with  standard” 


r 


Table  3.7 


Operating  Characteristics  of  a 72-Item  Inventory  System  Using 
the  Statistical  Power  Approximation:  Sensitivity  to  Input  Parameters 
(26-perlod  Revision  Interval,  26-perlod  Revision  History) 


INPUT  PARAMETERS 

Mean 

Tr/h  K/h  Lead  time  Demand 


OPERATING  CHARACTERISTICS 
E(COST) 

E (PERIOD-END  INVENTORY) 

E (BACKLOG  COST) 

E (BACKLOG  FREQUENCY) 

PROPORTION  OF  DEMAND  BACKLOGGED 
E (REPLENISHMENT  COST) 

E (REPLENISHMENT  FREQUENCY) 


■f  1 1 ill 
1 ^ 1 1 
+ + ? + + f 

+ + + 1 + + 

'I'  + •I'  + + 1 

1 + + 

+ + + + + 


+ + + + 1 + 

f + + 1 + + 

ill  111 

111  ? ? - 

111  ? 1 1 

111  111 

111  111 


POLICY  PARAMETERS 
D 
s 
S 


Ill 

111  111 

111  111 


111  111 
111  111 
111  111 


The  three  symbols  In  each  column  refer  to  Negative  Binomial  demand 
distributions  with  Standard  Deviation/Mean  - 1,  Variance/Mean  = 9,  and 
Variance/Mean  ■ 3,  respectively. 

1 indicates  a monotonlcally  increasing  relationship; 

1 indicates  a monotonlcally  decreasing  relationship; 

- indicates  an  approximately  constant  relationship; 

? indicates  that  no  simple  relationship  could  be  discerned. 


deviation- to-mean  ratio  of  1.  In  fact,  the  values  for  the  o/u  “ 1 
system  are  seen  to  be  quite  close  to  the  observed  values  for  the 
a^/p  “ 9 or  o^/p  ■ 3 systems. 
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Table  3.9 

Backlog  Frequencies  for  a Multi-Item  Systems 
with  Negative  Binomial  Demands 

I IT  / h 


9 

99 

1/d-Hi) 

.2 

.1 

.01 

Standard  Deviation/Mean  « 1 

.211 

.106 

.015 

Variance/Mean  ■ 9 

.205 

.107 

.019 

Variance/Mean  * 3 

.227 

.106 

.012 

I 
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The  manager  of  an  inventory  system  operating  with  incomplete 
demand  information  will  need  to  forecast  system  behavior,  especially 

to  justify  the  installation  of  scientific  control.  Forecasting  may 

I 

also  be  a routine  requirement  because  of  the  periodic  revision  of  the 
control  parameters  or  even  as  part  of  a regular  budgeting  procedure. 

Forecasting  methods  have  been  considered  in  MacCormick  (1974), 

Estey  and  Kaufman  (1975),  Ehrhardt  (1976),  and  Kaufman  and  Klincewicz 
(1976).  We  gather  forecasting  information  using  the  same  techniques; 
these  involve  double  use  of  the  demand  history,  for  control  and  fore- 
casting. Sample  distributions  of  the  forecasts  are  compared  with  those 
of  the  operating  characteristics,  which  were  obtained  in  Section  3. 

1 

The  major  Issues  to  be  resolved  here  are  the  extent  of  bias  and  the 
> level  of  dispersion  of  the  forecasts,  together  with  the  variation  of 

these  with  system  settings,  including  the  demand  process,  unit  costs, 
the  replenishment  leadtime,  the  revision  interval,  and  the  lengths  of 
the  histories  used  for  revision  and  forecasting.  Little  consideration 
is  given  to  an  Important  sampling  problem,  namely,  that  since  it  is 
usually  prohibitively  expensive  to  gather  a history  of  demands  for  every 
I item  in  a large  multi-item  system,  a sample  of  the  items  must  be  selected 

I to  make  a forecast.  The  accuracy  of  such  a forecast  is  likely  to  be 

i critically  dependent  on  the  number  of  items  in  the  sample  and  its  compo- 

sition. Forecasts  made  for  the  multi-item  system  in  this  simulation 
experiment  have  used  every  item  in  the  system.  Our  work  examines  the 
case  of  a revision  Interval  of  26  weeks  and  a revision  history  length 


of  26  weeks. 
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In  Che  simulation  experiments,  point  forecasts  of  the  expected 
values  of  system  operating  characteristics  are  made  each  time  the 
control  parameters  (s,S)  are  revised.  The  demand  history  used  to  make 
a forecast  has  been  selected  to  contain  demands  for  the  same  number  of 
periods  as  the  revision  interval. 

Average  values  of  the  forecasts  of  expected  costs  per  period  for 
the  a/y  * 1 , a^/y  = 9,  and  a^/;i  » 3 multi-item  systems  are  shown  in 
Table  4.1,  the  sample  standard  deviations  being  shown  In  parentheses. 

The  sample  statistics  for  forecasts  shown  In  Table  4.1  can  be 
compared  with  those  In  Table  3.1  for  average  costs  per  period.  The 
comparison  with  average  costs  per  period  is  explicit  in  Table  4.2,  which 
contains  the  percentage  differences  between  the  forecasts  and  the 
estimated  costs. 


Table  4.1 

Average  Forecasted  Values  for  Expected  Costs  per  Period  for  a Multi  Item 
Negative  Binomial  System  Controlled  with  Statistical  Power  Approximation 


Cost  Components 

Standard 

Deviation/Mean  = 1 

Variance/Mean  = 9 

Variance/Mean  = 3 

Inventory 

1959  (79) 

[65.5] 

1966(68) 

[65.6] 

1283(28) 

[56.2] 

Backlog 

384  (154) 

[12.8] 

405(57) 

[13.5] 

273(30) 

[12.0] 

Replenishment 

647  (14) 

[21.6) 

626(15) 

[20.9] 

727(11) 

[31.8] 

Total 

2991  (104) 

[100.0] 

2996(97) 

[100.0] 

2283(44) 

[100.0] 

Note:  Estimates  of  the  standard  deviations  for  the  forecasts  are  shown  in 
parentheses  and  each  cost  component's  percent  of  the  total  Is  shown  in 
brackets. 


ll 
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Table  4.2 

Average  Costs  Per  Period  for  Multi  Item  Negative  Binomial  System 
Controlled  by  Statistical  Power  Approximation:  Comparison  of 
Forecasts  with  Corresponding  Estimated  Expected  Values 


Cost  Components 

Standard 

Deviation/Mean  = 1 

Variance/Mean  = 9 

Variance/Mean  = 3 

Inventory 

1.5 

1.2 

1.1 

Backlog 

55.0 

59.3 

37.7 

Replenishment 

0.8 

1.2 

0.6 

Total 

14.5 

17.1 

7.5 

Note:  Table  4.2  shows  the  percentage  value  of  the  underestimate  by  the 


forecast  of  the  realized  value. 

On  the  average,  the  forecasts  underestimate  every  component  of 
expected  cost  for  all  systems,  a result  which  is  consistent  with  all 
previous  analysis.  The  underestimation  for  the  constant  standard  deviation/ 
mean  system  is  14.5%  which  compares  with  17.1%  for  the  variance/mean  =9 
system  and  7.5%  for  the  variance/mean  = 3 system.  Like  these  systems,  the  ; 

major  source  of  underestimation  is  the  backlog  cost  component  which  is  I ^ 

underestimated  by  55.0%.  We  note,  however,  in  Table  4.1,  that  for  the  standard 
deviation/mean  “ 1 system,  the  standard  deviation  of  the  backlog  cost  estimate  is 

1 : 

154,  which  Is  much  greater  than  the  values  for  the  variance/mean  “ 9 system  (57) 

i . 

or  the  variance/mean  ■ 3 system  (30).  | 

H 

Table  4.3  shows  the  percentage  difference  (bias)  between  the  fore- 
casted cost  components  and  the  observed  simulation  values.  For  the  j | 

constant  standard  devlatlon/mean  system,  the  bias  in  the  forecast  of  j j 

i 

expected  average  total  cost  and  the  backlog  cost  components  tends  to  j 


A 
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Increase  with  Increasing  mean  demand  and  the  bias  In  the  forecast  of  the 
replenishment  cost  component  tends  to  Increase  with  Increasing  leadtime. 
For  the  constant  variance- to-mean  systems  these  relationships  were  all 
monotonlcally  decreasing.  Other  trends  In  the  bias  of  cost  component 
forecasts  with  respect  to  Input  parameters  follow  the  same  pattern  In 
all  three  systems. 

Table  4.4  demonstrates  biases  for  forecasts  of  operating  character- 
istics other  than  costs,  expressed  as  percentages  of  the  estimates  of  the 
actual  expected  values  of  the  operating  characteristics.  We  see  that 
whereas  for  the  standard  devlatlon/mean  1 system,  bias  In  the  forecast 
of  backlog  frequency  and  proportion  demand  backlogged  tends  to  Increase 
with  Increasing  mean  demand  and  the  bias  In  the  forecasts  of  replenishment 
frequency  tends  to  Increase  with  leadtime,  these  same  relationships  are 
monotonlcally  decreasing  for  the  constant  varlance-to-mean  ratio  systems. 
The  only  characteristic  where  the  bias  In  the  forecast  Is  consistently 
worse  for  the  constant  standard-deviation-to  mean  system  is  period-end 
inventory,  but  overall,  the  bias  for  this  characteristic  Is  only  1.5% 
below  actual  and  should  not  Indicate  that, use  of  the  Power  Approximation 
algorithm  and  retrospective  simulation  for  a constant  standard-deviation- 
to-mean  system  results  In  any  significant  problems  which  are  not  also 
present  for  constant  varlance-to-mean  systems. 


Forecasting  for  a 72-Item  System  Under  Statistical  Control; 

Estimated  Percentage  Difference  Between  Forecasts  and  Expected  Average  Costs  per  Period 
[Revision  Interval  26  Periods,  Revision  History  26  Periods,  Forecasting  History  26  Periods 
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